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How to hide large scale outflows: size constraints on the 
jets of Sgr A* 



o 
o 

(N 
cn 



> 

cn 
\o 

(N 
O 
l> 

o 
o 



Sera Markoff^*, Geoffrey C. Bower^ and Heino Falcke^ 

^ Astronomical Institute "Anton Pannekoek", University of Amsterdam, Kruislaan 403, 1098SJ Amsterdam, the Netherlands 
^601 Campbell Hall, Astronomy Department & Radio Astronomy Lab, Berkeley, CA 94720, USA 

^ Department of Astrophysics, Radboud University, Postbus 9010, 6500GL Nijgmegen; ASTRON, Postbus 2, 7990AA Dwingeloo, the Netherlands 



ABSTRACT 

Despite significant strides made towards understanding accretion, outflow, and emis- 
sion processes in the Galactic Center supermassive black hole Sagittarius A*, the 
presence of jets has neither been rejected nor proven. We investigate here whether the 
combined spectral and morphological properties of the source at radio through near 
infrared wavelengths are consistent with the predictions for inhomogeneous jets. In 
particular, we construct images of jets at a wavelength of 7mm based on models that 
are consistent with the spectrum of Sgr A*. We then compare these models through 
closure quantities with data obtained from the Very Long Baseline Array at 7mm. We 
find that the best-fit jet models give comparable or better fits than best-fit Gaussian 
models for the intrinsic source found in previous analyses. The best fitting jet models 
are bipolar, are highly inclined to the line of sight {9 > 75°), may favor a position 
angle on the sky of 105°, and have compact bases with sizes of a few gravitational 
radii. 
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1 INTRODUCTION 



Sgr A* is the compact radio source in our Galactic center, 
origin ally discovered over 30 years ago by iBalick fc BrowrJ 
1 19741 ). Years of dedicated observations of stellar orbits 



(e.g., iGhez et all I2OO0I : ISchodel et all l2003f) and precise 
high-resolution radio astrometrv jBacker fc Srameld 1 19991 : 
iReid et al.|[2003t ) have led to Sgr A* being unambiguously 
associated with the central supermassive black hole. Most re- 
cently, the first mea surements of t he ins t rinsic size of Sgr A * 
have been achieved iBower et all l|2004l '): IShen et all l|2005h : 
iBower et al] l|2006h . giving information about spatial struc- 
tures extremely close to the black hole. 

For many years Sgr A* was only known to emit in the 
radio bands, with a flat/inverted spectrum fairly typical of 
the compact cores of other near by low-luminosity galaxies 
(e.g. IHoII1999| : iNagar et alll2002f ). However, the absence of 
infrared and higher energy emission was puzzling given that 
at least some nonthermal accretion activity would be ex- 
pected for a source that is powered by weak accretion. The 
first positive identification of Sgr A* in the X-ray band 
with the Chandra Obs ervatory did not immediately settle 
the issue (jBaganoff et al... 2001. . .2003, ). The dominant quies- 
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cent emission turned out to be extended and nonvariable, 
and thus likely associated with hot gas within the Bondi 
capture radius of the black hole (jQuata crt 2002.) . In con- 
trast, the approximately daily fiares of nonthermal X-ray 
emission discovered later seem to originate within tens of 
Tg = GM 1 1? from the black hole itself. Since this discovery, 
Sgr A* has also been identified in the near infrared (NIR), 
where it shows correl ated variability with the X-ray band 
on similar timesca les (|Genzel et al.ll200^ : iGhez et ai.ll2004l : 
lEckart et al.ll2003 ). While this suggests a low-level of active 
galactic nuclei (AGN)-like behavior, the luminosity of Sgr 
A* (~ 10~^LEdd) is weak enough to raise questions about 
comparisons with more luminous accreting black holes. 

Several models have been developed over the years 
to explain the broadband e mission of Sg r A*, rang- 
ing from Bondi-Hoyle infall l|Melial Il992f ). to various 



radia t ively-ineffici ent accretion fio ws fRIAFs: | Naravan et al.l 



19981: lYuan et a l' 2003), to jets jF alckc fc Biermannlll995l : 



Falcke fc Markofi .2000: iMarkoff et al.ii200ll ). and combina- 
tions thereof I Yuan et al.l 20021 ). The persistence of such a 



wide range of models can be attributed to some extent to the 
lack of constraints on the nonthermal part of the X-ray spec- 
trum. Fitting the most compact "submm bump" region of 
the spectrum results in fairly similar internal parameters for 
all current models, and this "theoretical degeneracy" cannot 
easily be broken without better morphological information 
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from Very Large Baseline Interferometry (VLBI). Unfortu- 
nately with current sensitivity and resolution limits, most 
structure in the source is washed out by a strong scattering 
medium in the central Galactic regions (see, e.g. Bower et al., 
I2OO6I ). 

Recently, however, several new observational techniques 
have been developed which may help discern between var- 
ious models. For instance, the stringent limits placed on 
the accretion rate [M ^ 10 ~^ — lO'^Ma/yr) by measure- 
ment s of lin e ar polarization jAitken et al.ll200(]|:lBower et al.l 
I2OO3I . I2OO5I : iMarrone etall l2006bl : iMacauart et al.l I2OO6I ) 
have ruled out classical versions of the Bondi-Hoyle 
and Advection-Dominated Accretion Flow (ADAF) mod- 
els. Similarly, better determinations of the frequency- 
dependence of the electron scattering l aw in the Galac- 
tic center (GC) fB ower et all I2OO4I : IShen et al.1 l2005l : 
iBower et ahl uOOG) have resulted in new constraints on mod- 
els via their size- versus- frequency predictions. While the 
different groups have found the index of the size-versus- 
frequency relation to range from ~ 1 — 1.6, clearly any suc- 
cessful model must be stratified (optically thick and thus 
having a photosphere whose observable size varies with fre- 
quency) to achieve this. The determination of the scattering 
law to a high degree of accuracy has allowed, for the first 
time, a dependable measurement of the intrinsic size of Sgr 
A* along one axis as a function of frequency. This break- 
through, along with the expectation of eventually determin- 
ing the size in the other axis, means we are finally at a key 
point where differences between models can be empirically 
tested. 

In this paper, we use both the spectral data in combi- 
nation with the new VLBI measurements of th e source pho- 
tosphere at 43 GHz (from iBower et al.l |2004| plus one new 
observation, see below) in order to place new constraints 
on jet models. In Section [2j we expand on the motivations 
for this project, in Section [3] we introduce the model, in 
Section U we explain the methodology and summarize our 
results in Section[5l and discuss our conclusions in Section[6l 



2 THE EVIDENCE FOR JETS IN SGR A* 

Because no jet in Sgr A* has yet been directly imaged, it is 
important to first discuss the evidence in favor of jets in Sgr 
A*. The lack of a resolved core/jet structure is not surprising 
given the low luminosity of Sgr A*, which suggests a small 
angular size for the jet, and the scattering screen in our line 
of sight towards Sgr A*, which obscures small structures. 
Previous modeling of the structure of Sgr A* has succeeded 
in separating the intrinsic and scatter-broadened images of 
Sgr A* via a Gaussian parameterization of the intrinsic size. 
A primary goal of this paper is to go beyond this simple 
parameterization. 

In fact, there are several strong arguments for jets in 
Sgr A*. On a purely theoretical level, some form of jet pro- 
duction seems to go hand-in-hand with accretion around 
black holes, both at the galactic as well as stellar scales. 
In stellar black holes accreting from a binary companion, 
or X-ray binaries (XRBs), jet production is observed to 
be cyclic over outburst cycles. The strongest (relative to 
the system energetics) and steadiest jets occur during the 
low-luminosity state, called the Low/Hard State, while dur- 



ing the highest lum inosity state, the jets appear quenched 
dFender et al.lll999l ). The low-luminosity jets are compact 
and self-absorbed with a flat/inverted spectrum, and corre- 
lated radio/X-ray variability has demonstrated that the jets 
increasingl y dominate the po wer output as the luminosity 
decreases ( Fender et al.|[2003 ). The weakest accreting black 
hole we can study with reasonable statistics besides Sgr A* 
is the XRB A0620 -00, in which r adio emission has recently 
been detected (jGallo et al.ll2006l ). At an X-ray luminosity 
of < 5 X 10~^ I/Edd, very close to that of Sgr A*, efficient 
jets are still produced in this black hole, with characteris- 
tics matching those at higher powers. If general relativity's 
basic prediction of scaling black hole physics holds, this is a 
strong argument for jet production in Sgr A*. 

The radio spectrum, radio variability, and high- 
frequency linear polarization are all similar to other nearby 
low-luminosity AGN (LLAGN; Ho 1999; Falcke fc Bie rmann| 
1 19991 : IBower et al.ll2002l : iNagar et al.ir20oi l2005l )71vlost of 
the observed cores are accretion-powered, and have the sig- 
nature flat /inverted, sel f-absorbed radio spectrum associ- 
ated with compact jets (|Blandford fc Konig]||l979l ). While 
the jets can only generally be resolved in the brightest 
sources, when they are resolved they dominate the unre- 
solved core by at least a factor of a few. The results of these 
surveys strengthen the arguments for a jet in Sgr A* based 
on its radio spectrum and polarization. 

One source that is particularly interesting because of 
its many parallels with Sgr A* is the nucleus of the nearby 
LLAGN M81. M81* is our nearest LLAGN besides Sgr A*, 
and resides in the same kind of spiral galaxy as the Milky 
Way. Its mass has been derived from line spectroscopy (using 
HST; Devereux et al. 2003) to be 7 x lO'' Mq, only ~ 30 
times the mass of Sgr A*. M81* also possesses th e typical 
compact fiat/inverted core spectrum (|Falckdll99d : Markoff 
et al., in prep.) and, more importantly, the same high levels 
of circular rather than linear polarization in the cen timeter 
radio band as Sgr A* (|Brunthaler et al.ll200ll . I2OO6I '). 

The M81* jet is one-sided, very small (700-3600 AU 
depending on the frequency, with a roughly ~ l/i/ de- 
pen dence), and exhibits oc casional bends in its morphol- 
ogy (iBietenholz et allbOOOl ). Scaling the size by mass alone 
would argue for a ~ 20 — 120 AU jet in Sgr A*, but the 
observed size should also scale with luminosity, depending 
on the particulars of the jet model and frequency. In fact, 
the jet nature of M81* was difficult to establish due to the 
high level of compactness. Taking into account Sgr A*'s five 
orders of magnitude lower power, as well as the scattering 
screen, it is not surprising that no jet has yet been detected 
in our Galactic center. 

Another argument in favor of jets comes from the recent 
detection of short time delays of about 0.5-1 hr between 43 
and 22 GHz for waves of va riability trave ling from high to 
low frequencies tYuscf-Zadc h et al] |2006| ). This variability 
is fully consistent with outfiowing, adiabatically expanding 
blobs of plasma, as would be expected for jets (in fact, the 
model the authors use to interpret their results was devel- 
oped in this context). 

Finally, the recent size-versus-frequency scaling detec- 
tions support an optically thick, stratified model such as 
a self-absorb ed jet. While the predict ions of the jet model 
presente d in iFalcke fc Markofj l|2000l ), as well as that of 
RIAFs l|Yuan et aLlbOOd ). are consistent with a l/i/ seal- 
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ing, and thus wi th the results in iBower et alj l|2004 ) and 
IShen et ahl l|2005l ). they di sagree with the stee per index de- 
termined more recently by lBower et al.l (|2006l ). If this latter 
result is indeed correct, it suggests that the current versions 
of all models, jets included, need to be modified to show a 
stronger dependence on observing frequency. Because this 
issue is still under debate, however, in this paper we are still 
using the original scaling relation. 

Although the circumstantial evidence is significant, 
there are other complications which could argue against jets. 
For instance, XRBs in their steady-jet producing Low/Hard 
state display a correlation between their radio and X-ray 
luminosities that hol ds over at least seven orders of mag - 
nitude in luminosity (jCorbel et al.ll2003l : [Callo et al.ll2003h . 
Among other things, this correlation can be used as a gauge 
for "typical" levels of activity. The recent radio detection 
of A0620-00 falls exactly on the correlation, extending it to 
even lower luminosities and indicating that the same mech- 
anism is at work as in brighter sources where jets can be 
imaged. If the physics driving the correlation scales in a pre- 
dictable way with mass, it should apply to LLAGN as well, 
where the mass enters mainly as a normalization factor for 
the same correlation slope. This relationship between radio 
and X-ray luminosities and mass is called "the fundamental 

plane of black hole accretion" and has been explored in sev- 

I 1 1 ) I ' 1 1 — — — I 

eral recent pa pers (Merloni et al. 2003; Falck e et al.l 120041 : 

iKording etal.. .2006; Merloni e t al. 2006). When Sgr A* in 

quiescence is placed on this plane, it falls well below the 

correlation in predicted X-rays, given its radio luminosity. 

One could interpret this as the complete dominance of the 

jet over infiow processes at the lowest of luminosities, but it 

could also mean that the emission mechanisms themselves 

have undergone a transition to a different mode of emission 

entirely. 

In order to try to cast new light on these long-standing 
ambiguities and place more stringent constraints on the pos- 
sible presence of jets in Sgr A*, we have developed a new 
method to combine spectral and morphological data. Our re- 
sults will set the stage for future tests with upcoming VLBI 
observations at millimeter and submillimeter wavelengths, 
where the morphology is less affected by scattering and res- 
olution is comparable to a few rg. 



3 MODEL 

Like most models involving optically t hick, coUimated out- 
fiows, we build on the initial work of iBlandford fc Konigll 
l| 19791 ). These authors demonstrated the "conspiracy" of 
how a perfectly fiat spectrum {a ~ 0, oc can result 

from a superposition of self-absorbed contributions along a 
conical, idealized jet. When more realistic physics such as 
bulk acceleration, full particle distributions and cooling are 
included, compact jets show a slight spectral inversion in the 
radio wavebands, with a ~ 0.0 — 0.2. The model used here is 
based on a model developed for Sgr A* (jPalcke fc Marko3 
l2000h . which has been significantly modified to extend to 
XRBs and LLAG N in general. For a d etailed description 
see the appendix in lMarkoff et al.l (|2005l ): we provide only a 
brief summary below. 

The model is based upon four assumptions: 1) the to- 
tal power in the jets scales with the total accretion power at 



the innermost part of the accretion disk, Mc^, 2) the jets are 
freely expanding and only weakly accelerated via their own 
internal pressure gradients, 3) the jets contain cold protons 
which carry most of the kinetic energy while leptons domi- 
nate the radiation and 4) particles have the opportunity to 
be accelerated into power-law tails. In sources accreting at 
higher levels this latter point would be more important, but 
as we will show later, there is not much capacity in the Sgr 
A* spectrum for significant particle acceleration. 

The base of the jets consist of a small nozzle of constant 
radius where no bulk acceleration occurs. The nozzle absorbs 
our uncertainties about the exact nature of the relationship 
between the accretion flow and the jets, and fixes the initial 
value of most parameters. Beyond the nozzle the jet expands 
laterally with its initial proper sound speed for a relativistic 
electron/proton plasma, 7s/3sC ~ 0.4c. The plasma is weakly 
accelerated by the resulting longitudinal pressure gradient 
force, allowing an exact soluti on for the ve locity profile via 
the Euler equation (see, e.g., lFalckeiri99d ). This results in 
a roughly logarithmic dependence of velocity upon distance 
from the nozzle, z. The velocity eventually saturates at large 
distances at Lorentz factors of Fj >2-3. The size of the base 
of the jet, ro, is a free parameter (but expected to fall within 
several rg) and once fixed determines the radius as a function 
of distance along the jet, r{z). There is no radial dependence 
in this model. 

The model is most sensitive to the fitted parameter Nj , 
which acts as a normalization. It dictates the power ini- 
tially divided between the particles and magnetic field at 
the base of the jet, and is expressed in terms of a fraction 
of the Eddington luminosity L-^dd = 1.25 X lO^^Mbh srg 
s~^. Once TVj and ro are specified and conservation is as- 
sumed, the macroscopic physical parameters along the jet 
are determined. We assume that the jet power is roughly 
shared between the internal and external pressures. The ra- 
diating particles enter the base of the jet where the bulk 
velocities are lowest, with a quasi-thermal distribution. In 
higher power jets, a significant fraction of the particles are 
accelerated into a power-law tail, however in Sgr A* this 
seems to be less of an effect. The particles in the jet ra- 
diatively cool via adiabatic expansion, the synchrotron pro- 
cess, and inverse Compton upscattering; however, adiabatic 
expansion is assumed to dominate the observed effects of 
cooling. Because S gr A* has no "s tandard thin accretion 
disk" (e.g. .Shakura fc Sunvaevlll973l ). nor even a fossil disk, 
which would be apparent in the infrared l|Falcke fc Melial 
|l993), the photon field for inverse Compton upscattering 
is entirely dominated by locally produced synchrotron pho- 
tons. Fig. [T] shows an example of the resulting broadband 
spectrum plotted against the data for Sgr A*. 

Besides those mentioned above, the other main fitted 
parameters are the ratio of length of the nozzle to its ra- 
dius ho, the electron temperature Te, the inclination angle 
between the jet axis and line of sight 6i and the equiparti- 
tion parameter between the magnetic field and the radiating 
(lepton) particle energy densities, k. 

Aside from Sgr A*, this class of mo del has been suc- 
cessfully applied to several LHS XRBs (.Markoff et al.ll200ll . 
I2OO3I : iMarkoff et al.ll2005l : iMigliari et all" 20071: GaUo et al., 
in prep.) and other LLAGN (|Yuan et al.ll2002l : Filho et al., 
in prep., Markoff et al., in prep.). As would be expected from 
the existence of the fundamental plane, all significantly sub- 
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log,„ u (GHz) 

Figure 1. The broadband data set for Sgr A* used to constrain 
our models , taken from tlie aver age spectrum up until tiie sub- 
millimeter jMelia fc F al ck e)|2p0ll). with additional low frequency 
points from iNord et al.1 ll20o3rand Roy & Pramcsh Rao (2004) 
and i nfrared data from ICenzel et al.l ||2003|) and iGhez et al.l 
l|2004l ) . The X-ray "bow-ties" represent the quiescent (lowest), 
average daily Chandra flare (midd le) and brightest Chandra (top) 
powe r-laws with errors indicated ijBaganoff et al.ll200ll : iBaganog 
I2OO3I) . Th e "V" shape indicates the two XMM-Newton flares pre- 
sented in iBelanger et al.l I^OO^) • The solid curve shows a repre- 
sentative quiescent model with synchrotron and synchrotron self- 
Compton peaks. The dotted lines illustrate the contribution of the 
quasi-thermal particles from each increment along the jet, which 
superimpose to give the characteristic flat/inverted synchrotron 
spectrum. 



Eddington accreting black holes do seem to share some basic 
underlying physics across the mass scale. However, as men- 
tioned above, Sgr A* does not participate in the radio/X-ray 
correlation and can only be reconciled into this picture if sig- 
nificant particle acceleration is lacking. This is a very inter- 
esting point, because the appearance (or non-appearance) 
of a jet is strongly dependent on its internal particle distri- 
butions. A power-law tail of accelerated particles results in 
more optically thin synchrotron emission over a broader fre- 
quency range from each jet increment. Thus when observing 
at a single frequency, a larger range of increments are able 
to contribute to the profile, resulting in a larger jet image, 
Eis we show in Fig. O 



4 METHODOLOGY 

4.1 Modeling the spectral data of Sgr A* 

In order to break the current degeneracy based on modeling 
the spectral data alone, we are for the first time calculating 
the predicted size and morphology of jet models which give 
a good description of the broadband data. After convolving 
these "maps" with the scattering screen, we then compare 
the results to closure quantities from VLBI, which provide 
information about the structure. The model predictions are 
analyzed in the same way as the observations. 

The 41 models presented here were chosen to represent 
a range which samples the full parameter space within the 
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Figure 2. Illustration of the effect of particle acceleration on 
the observed jet profile. Panel a) shows the quiescent model 3, 
which provides the best statistical description of the radio data. 
Panel b) shows a model with the exact same parameters except 
for the addition of significant particle acceleration, where 75% 
of the particles in the quasi-thermal distribution are accelerated 
into a power-law with energy index 2.2, and a cutoff 50 times 
higher than the minimum energy (which is fixed to the peak of 
the thermal distribution). Panels c) and d) show the profile this 
model produces on the sky, prior to convolution with the scatter- 
ing Gaussian. The scale of the images is 1 mas. The images are 
saturated to emphasize the brightest (darkest) parts of the jets. 
A tail of accelerated particles serves to significantly lengthen the 
jets profile at a given frequency. 



constraints of a reasonably good (x^ < 1) description of the 
radio through IR. For quiescent models, they are further 
constrained to fall within a factor of a few under the X- 
ray quiescent limits from Chandra. Flaring models are those 
which can account for either the average daily flaring flux 
or the highest detected with Chandra in the X-ray band, 
via some form of heating, accelerating or increased power 
compared to the quiescent state. We initially conducted a 
very rough search of a large region of parameter space for 
the first 20 models, then focused on a smaller region to ex- 
plore properties nearest the best-fitting models, as well as 
for flares. A systematic exploration of the parameter space 
would be at least a 6 dimensional parameter cube, which 
we deemed too computationally intensive for this initial test 
study. 

Because there is still significant uncertainty about what 
to consider the "quiescent" versus "flared" IR flux amplitude 
and slope, we did not include the IR in the calculation 
but rather required the model to fall reasonably within range 
of the errors. The inability to constrain the exact IR and 
X-ray quiescent flux accounts for almost all of the allowed 
range in the fitted parameters for the quiescent state, other- 
wise the parameters would be fairly tightly determined. In 
this way the addition of morphological fitting can constrain 
the quiescent contribution above the submm band. We also 
compared our models to the simultane ous radio through IR 
data set presented in lAn et all l|2005h . and found the level 
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agreement between the two data sets so good that we did 
not feel it provided an additional constraint. One of our sec- 
ondary goals was to explore observable differences in the 
photosphere during flaring states. 

Table [1] lists the models and their parameters, along 
with relevant comments. 



Table 1. Jet model parameters 
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PL:p = 3.4, 7o,max = 2 X 10^ 
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Q 
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y 
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34 


F 


19 
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1.0 


85 


1.0 


15 


y 


2 


Zacc = 10, p = 1.7, u/f = 0.014, plf= 1 X 10^" 
Biggest Chandra flare, synch. 


10.35/5 


35 


F 


11.5 


2.5 


0.95 


85 


1.0 


50 


y 


2 


Zacc = 10, p = 1.6, u/f = 0.014, plf= 6 X 10^" 
Av. Chandra flare, synch. 


5.90/5 


36 


F 


50 
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85 
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y 


2 


PL:p = 2.3, To, max = 500 
Biggest Chandra flare, SSC 
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PL:p = 1.01, To, max = 220 
Zacc = 5, p = 1.01, u/f = 3 X 10"'', plf= 5 X 10"* 
Av. Chandra flare, SSC 
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Av. Chandra flare, SSC 
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Biggest Chandra flare, SSC 
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41 


F 


25 


2.5 


1.01 


85 


2.0 
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y 


2 


Zacc = 5, p = 2.7, u/f = 3 X 10~^, plf= 0.04 


4.193/5 



Steep XMM fiarc, synch. 



Note. — This column describes other adaptations to the standard model. "/3 stretched" means that we increased the velocity as a function of distance along 
the jet by a factor depending on that distance. The other comments refer to various ways of accelerating particles in the jets. For rows with four additional 
parameters, Zacc is the location of the acceleration region, p is the particle index, u/f are plasma param eters which determi ne the rate of acceleration, and 
plf is the fraction of particles accelerated out of the original quasi-thermal distribution (see Appendix in iMarkoff et al.||2005l for details). For rows with two 
additional parameters, the particles are assumed to be accelerated already in the nozzle, in a power law with 7^ rnin corresponding to the input temperature, 
particle index p and maximum lepton Lorentz factor TTiaa:- 

The statistic is sensible for quiescent models only, since it is calculated using an averaged quiescent spectrum. We include the value for flare models just 
as a reference. 
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4.2 Analysis Technique 

The jet emission is calculated along its length in increments. 
In order to determine the appearance of the jet on the sky, 
we calculate the contribution to the synchrotron spectrum 
at 43 GHz from each increment, assumed to be evenly dis- 
tributed over the radi us and increment width. Relativistic 
angle aberration (e.g. iLind fc BlandfordI [l98^ for the in- 
crements' bulk Lorentz factors is taken into account. This 
"profile" is then fed into an IDL routine which creates a 
FITS image of the jet. Each model was then rotated by po- 
sition angles in steps of 15 degrees covering the full range of 
angle. Furthermore, once the jet was placed with the speci- 
fied rotation on the image, we convolved it with a Gaussian 
ellipse of the scattering as determined below. 

Jet models were imaged on a 2001 x 2001 grid with a 
pixel resolution of 14 /^arcsec. Fig.[3]shows the underlying jet 
model and the resulting scatter-broadened model in linear 
and logarithmic scales. The large-scale differences seen in the 
logarithmic representation do not make a significant contri- 
bution to our ability to differentiate between these models, 
since the total flux density in the outer regions is very small. 

We directly compare the jet models with high resolution 
data obtained at a wavelength of 43 GHz (7 mm). The data 
are obtained primarily with the Very Long Baseline Array 
(VLBA) and in some cases include a single Very Large Array 
( VLA) antenna . Eigh t epochs of observations are described 
in lBower et al.l (|2004h . In addition, we include new observa- 
tions obtained with the VLBA and the 100m Green Bank 
Telescope on 18 May 2004 (experiment code BB183). These 
observations were reduced in the same method as the earlier 
epochs with calibration for single-band delay and multi-band 
delay and rate. 

We construct closure amplitude and closure phase from 
the visibility data. The closure phase is the sum of interfer- 
ometric phases for a triangle of baselines. The closure am- 
plitude is a product of interferometric amplitudes for base- 
line quadrilaterals. Analysis of the closure quantities is less 
sensitive than the analysis of calibrated visibilities because 
of the reduced number of degrees of freedom. The closure 
quantities are independent of amplitude and phase calibra- 
tion, however. This property which makes them valuable es- 
timators of source structure that are unbiased by systematic 
errors in calibration. 

In lBower et al.l l|2004l ') elliptical Gaussian models were 
fitted to the closure amplitudes for data sets at wave- 
lengths from 7 mm to 6 cm. This fitting produced a best- 
fit elliptical Gaussian as a function of wavelength. Com- 
bining the VLBI measurements with new measurements 
of the size at wavelengths between 17 and 24 cm based 
on VL A observations, a si ze-wavelength relation was deter- 
mined l|Bower et al.| [200^ . The scattering ellipse from the 
long wavelength observations was computed to be 1.31 x 0.64 
meis cm~^ in position angle 78°. The size of the ellipse scales 
as the wavelength-squared. Deviations from the wavelength- 
squared law at short wavelengths are indicative of the in- 
trinsic size becoming comparable to the scattering size. The 
magnitude of the scattering ellipse is determined by the 
spectrum of turbulent electron density fiuctuations. The ori- 
entation and axial ratio of the scattering ellipse are deter- 
mined by the magnetic field properties of the plasma in 
which the scattering originates. 



6 



30 



34 



41 




Figure 3. Images of models 6, 30, 34, and 41. In the left col- 
umn, we show the jet model with a linear stretch. In the middle 
column, we show the jet model after it has been convolved with 
the elliptical Gaussian representing scattering, also with a linear 
stretch. In the right column, we show the convolved jet model in 
a logarithmic stretch. The scale bars in the top row indicate 1 
milli-arcsecond. 

Each model was directly compared with closure quan- 
tities from the data. First, we added a noise bias to each 
model image equal to the best-fit value determined from 
elliptical Gaussian fitting for the data set. We also experi- 
mented with using a range of noise biases that went from 
to 2 times the best-fit value. We found that the minimum 
from this procedure was comparable to the for the 
best-fit noise bias. Second, we constructed the image two- 
dimensional FFT, which is the visibility plane representa- 
tion of the data. Third, closure quantities were computed 
for each model for the time and antenna sampling of the 
data set. Finally, reduced was computed for closure am- 
plitudes and closure phases for each model and each data 
set. 

In addition to jet models, we also created a model image 
for an elli ptical Gaussi a n tha t represents the best-fit Gaus- 
sian from iBower et ahl l|2004l ). The reduced xt = 1-9 from 
this fit is the baseline result that jet models must meet or 
surpass in order to remain viable. 

To demonstrate the ability of our method to discrimi- 
nate between models, we substituted the closure quantities 
from the data with closure quantities derived from model 41 
in three different position angles (90, 120, and 180 degrees). 
We then compared the substituted closure quantities with 
closure quantities from all models and position angles (Fig- 
ure |4]) . We computed the results for three different values 
of the noise bias. These results show that we can differenti- 
ate between position angles and models in the case of high 
signal-to-noise ratios (SNR). 

We also considered whether there are systematic differ- 
ences in the model between different data sets (Figure [S}. 
Seven of nine data sets are essentially consistent with each 
other. Data set BB130C shows a flat a function of 

position angle. This is consistent with larger than average 
noise (Figure |4)| , which was also seen in poor limits from 
the Gaussian fitting (Bower et al. 2004). Data set BS055C 
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Figure 4. Total xt ^ ^ function of position angle for models 6, 
30, 34, and 41 for the case where the data is replaced with closure 
quantities calculated from model 41 in position angle 90 deg (left 
column), 120 deg (middle column), and 180 deg (right column). 
Model closure quantities were computed for three different esti- 
mates of the noise, with the solid line representing the least noise, 
the dot-dashed line the middle case, and the dashed line the most 




1 50 200 250 

Position Angle (deg) 



Figure 5. Total xf, a-s a function of position angle for model 41 
showing each radio data set individually. 



shows a similar profile in Xi^ versus position angle but sig- 
nificantly larger values than average. This suggests that we 
may have underestimated the noise for this experiment. We 
have therefore dropped these two outlier experiments from 
all further modeling results. 



5 RESULTS 

In Figure [6] we show all of the model images used in the 
analysis. In Figure [7l|8l and [9] we show the closure phase, 
closure amplitude, and total Xv as a function of position 



Figure 6. Images of all models prior to convolution with the 
scattering ellipse, with a linear stretch. The scale for each image 
is 1 milliarcsecond. 
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Figure 7. Closure phase x5 as a function of position angle for all 
models. The dotted line represents the reduced the best-fit 

Gaussian model. 



angle for each of the models. In order to see details for the 
best-fitting models, we plot only on a scale of to 5. For 
several models, > 5i thus where no curve is present, the 
model is already strongly rejected. 

For a number of models, the minimum is less than or 
comparable to the best-fit Gaussian model. For all cases pre- 
sented here, Xv never achieves a significantly smaller value 
than the best-fit Gaussian model, which would allow un- 
equivocal rejection of that model in favor of a jet model. 
Instead, these results demonstrate that we can adequately 
but not uniquely model the data as a bipolar, relativistic 
jet. This result alone shows that jets in Sgr A* cannot be 
ruled out on the basis of their being unresolved. 

We see that most of the deviations in as a function of 
position angle are visible in the closure amplitude. For many 
models, the closure phase x^ is independent of position an- 
gle and is comparable to the value from Gaussian fitting. 
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Figure 8. Closure amplitude a function of position angle 

for all models. The dotted line represents the reduced for the 
best-fit Gaussian model. 
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Figure 9. Total as a function of position angle for all models. 
The dotted line represents the reduced for the best-fit Gaussian 
model. 

The closure phase results indicate that any deviations from 
axisymmetry in the source image are very small. The domi- 
nant role of the closure amplitude in variations with position 
angle is indicative of sensitivity to the size of the source in 
a given direction. 

In the East- West direction (90/270°), we have the best 
size determination because of the better resolution of the 
telescope in this direction. One can think of this as the 
data having smaller "error bars" around these angles. On 
the other hand, the scattering angle is largest, so asym- 
metries and extended components may be more obscured. 
In the North-South directions, the resolution is worse by a 
factor of ~ 3. Thus, minima in 0, 180 and 360 in 

several models are not significant. In these cases the over- 
all fit is very bad (as can be seen in the regions of better 
resolution around 90/270°) and the dips represent instead 



regions where our data are the least constraining. However 
in several models we see minima which are clearly offset 
from 90/270°, such as model 41 where the minimum occurs 
at ~ 105° (this can be most easily seen in Fig. [12}. The 
peaks seen at 0/180/360° suggest that even with the poor- 
est resolution, the asymmetry is too great to be consistent 
with those directions. The fact that the model is minimized 
at an angle where our ability to constrain the asymmetry 
is greater is suggestive, and its total x2 is in fact slightly 
lower than the best-fit Gaussian value. However, it is far 
too preliminary to claim a detection of a preferred position 
angle. These results do suggest, however, that with better 
resolution, especially in the N-S direction, the position an- 
gle of a jet may be constrained, particularly during fiaring 
states. Furthermore, many position angles are clearly ruled 
out, never achieving close to minimum for any spectrally 
consistent model. 

In Figs. I10H14I we show "scatter plots" of the mini- 
mum x^ from Fig. associated with some model parame- 
ters, for the 30 quiescent models only. The size and darkness 
of the circle/ellipse are inversely proportional to the x'^i i-e, 
large and black circles/ellipses are the best fits while lighter, 
smaller regions are not. The two data sets which are clearly 
discrepant from the others as discussed above, BB130C and 
BS055C, are not included. 

Fig. [To] demonstrates the clear selection of compact jets 
(whose smallest scale is the nozzle radius ro) and high in- 
clination angles. A much more stringent constraint than the 
spectrum alone is the combined effect of these two parame- 
ters on the jet profile. While any compact nozzle less than 
several rg, or inclination above ~ 75°, is acceptable spec- 
trally, the high level of symmetry required by the VLBI 
data strongly favors the most compact jets which are the 
most beamed out of the line of sight. Because the jets 
are mildly accelerating, the beaming-induced "dimming" in- 
creases along the jet axes, thus emphasizing the less elon- 
gated nozzle regions. These results are also a reassuring 
confirmation because it would be surprising and somewhat 
alarming if the jets were so misaligned as to be pointing sig- 
nificantly towards the Galactic plane in which we roughly 
sit. Fig [TT] also compares two geometrical parameters, this 
time the position angle on the sky versus the inclination an- 
gle. The best fit jet is therefore almost perpendicular to us, 
with a position angle on the sky of ~ 105° . Fig. [12] shows 
the clear peak in 1/x'^ at this angle. 

Fig. [I3]gives an example of how the additional morphol- 
ogy comparisons can also help constrain internal jet param- 
eters such as the equipartition of energy and electron tem- 
perature. While the overall range of spectrally-allowed tem- 
peratures spans a decade in temperature, the upper range 
clearly does not provide a compact enough jet profile. The 
equipartition parameter however is best constrained by the 
spectral fitting, which has already selected a rather narrow 
range. Values > 1 are magnetically dominated. 

Finally in Fig. [14] we show that some parameter degen- 
eracies clearly remain despite our new approach. Here we 
plot the electron temperature against the jet normalization 
parameter Nj. A clear range of acceptable values exists in 
both parameters, demonstrating for instance how a higher 
temperature electron distribution can compensate for lower 
power because of its more energetic emission. This can be 
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Figure 10. For each of the 30 quiescent models, we plot the 
minimum in PA from Fig. |9] for the indicated nozzle radius ro 
and inclination angle. The radius and greyscale (from white to 
black) are linear in l/x'^i f^nd smaller (larger circle) fits are 
plotted last. The largest, darkest circles have Xmin 
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Figure 11. Color scale the same as Fig. 1101 with symbols now el- 
lipses (axes scaled linearly inl/x^)to better illustrate the param- 
eter space, showing ^ ^ function of model PA and inclination 
angle. 



understood from the critical synchrotron frequency relation- 
ship Us oc Bj^, where oc Nj and 7^ oc Te. 



6 DISCUSSION AND CONCLUSIONS 

The most important conclusion of this paper is that a jet 
model, with reasonable physical assumptions about its ge- 
ometry and internal physics, is capable of explaining the ra- 
dio through IR (and higher, during flares) spectrum of Sgr 
A* and not be visible at all to us as an object with jet-like 
morphology! Aside from the overall low jet power, the lack 
of significant particle acceleration implied by Sgr A*'s IR 
spectrum would predict extremely compact jet profiles. Our 
results demonstrate that the lack of an imaged jet in Sgr 
A* is not necessarily due to any absence, but rather stems 
from a very compact, weak source combined with the rather 
extreme scatter broadening by Galactic electrons, and limits 
on our spatial resolution, especially in the N-S direction. 
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Figure 12. For all models, as a function of jet PA. The best 

fit occurs for 105/285° . 




5x10 



10" 1.5x10" 2x10" 2.5x10" 3xl0" 
Temperature (K) 



Figure 13. Same symbol definitions as Fig. 1101 showing x^ as 
a function of electron temperature and equipartition parameter 
(between magnetic and radiating particle energy densities, with 
> 1 meaning magnetically dominated). These are the two most 
important internal rather than geometrical parameters. 
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function of electron temperature and jet normalizing power (re- 
lated to, but slightl y less than, the total power; see Appendix in 
iMarkoff et n]|2005l') . There is a clear relationship between these 
two parameters. 
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However, even without being able to detect a fully elon- 
gated structure, the combination of spectral fitting with con- 
straints from comparison with VLBI morphology can signif- 
icantly limit the acceptable range of parameter space for jet 
models of Sgr A*. Figures [T0lH4l visually demonstrate these 
new limits, which are successful despite the preliminary ex- 
ploration of all parameter space. 

Not altogether surprisingly, the additional inclusion of 
size constraints from VLBI places tighter limits on the model 
geometry. For instance, while fitting the quiescent spectrum 
can only limit the inclination angle to > 45° , the addition of 
VLBI data in indicates a very clear preference for di ~ 90° . 
Similarly, the size constraints from VLBI also narrow the 
range in jet nozzle size from < 8rg to < 5rg with the best 
fits at the smaller end. In the context of jet models, this 
would require jet launching to occur very close to the black 
hole, within the innermost stable circular orbit (ISCO) for 
a Schwarzschild black hole. 

Size constraints also indicate that VLBI data can al- 
ready begin to constrain the orientation of jets on the sky. 
The best fit is found over a narrow range 90 — 120° centered 
at ~ 105°, in a region where the resolution of the VLBI is 
good enough to begin discerning the asymmetry. The pre- 
ferred PA is interesting, in that it could be related to the 
average position angle of the electromagnetic fields and thus 
give further clues about jet geom etry. Recent observa tions 
of variable lin e ar pola rization by iBower et al] Hooi) and 
iMarrone et all l|2006al ) observed PA changes of 30° - 60° 
over timescales of days to months. Infrared measurements 
of the polarization during flares also show significant vari- 
ability (|Trippe et al.l [20071 ). The variability is most likely 
intrinsic, although there may be a favored or mean intrinsic 
polarization PA in the various wavelengths, though currently 
they do not seem to agree with each other. Confirming both 
angles may ultimately provide important information about 
the helicity of the magnetic field threading the jets, or near 
the black hole. 

In addition. Figs. I7I9I clearly indicate a dramatic dif- 
ference in the goodncss-of-fit between quiescent and flaring 
models for Sgr A*. T his is because the mechanisms involved 
in creating the fla res (jMarkoff et al.ll200ll : lLiu fc Meha|[200^ : 
I Yuan et al.l 120031 ) are either heating or accelerating the ra- 
diating particles, which alters the optical depth and changes 
the jet profile on the sky. Our results strongly argue for fur- 
ther simultaneous X-ray and VLBI (eventually preferentially 
in the millimeter regime) monitoring of Sgr A*, where these 
methods can strongly limit the contributions of acceleration 
and heating, respectively. 

In conclusion, we find that the combination of broad- 
band spectral and morphological constraints gives encourag- 
ing and interesting limits on jet models (or any model) which 
cannot be obtained by spectral fitting alone. In particular, 
the current difficulty in constraining the high-energy contri- 
bution of the jets because of the dominant quiescent ther- 
mal X-ray emission highlights the need for new approaches. 
Including constraints from VLBI images offers a powerful 
method to break the current degeneracy in theoretical mod- 
els for Sgr A*'s emission, as well as better constraint indi- 
vidual models themselves. 

At 43 GHz and below, the key outstanding problem is 
to measure the two dimensional structure of Sgr A*. This 
requires a careful selection of North- South baselines that are 



sensitive to structure on the scale of a few hundred micro- 
arcseconds. However, it is important to note that electron 
scattering still acts to symmetrize the data at 43 GHz, thus 
mm/submm VLBI could be even more revealing for these 
types of studies. The advantage may, however, be offset by 
the fact that higher frequencies probe even smaller scales in 
the jets, which would be predicted to be as symmetric as an 
accretion flow. On the other hand, mm/submm VLBI will 
bring us to scales comparable to those probed by the IR/X- 
ray flares. While this will allow us to better observe simul- 
taneous flares in all three frequency bands, it also raises the 
question of how to distinguish the base region of a jet from 
an accretion flow. By identifying structural changes in mor- 
phology with spectral changes in a flare, the approach pre- 
sented in this paper will be able to constrain the geometry, 
particle distributions and emission mechanisms contributing 
to the flares. 

Finally we emphasize that the results presented here do 
not include modiflcations due to general relativistic effects 
near the black hole. For the current resolutions this may not 
be critical, but as we probe closer to the innermost regions 
with higher f requencies, this cl e arly needs to be taken i nto 
account (e.g. iFalcke et al.ll2000l : iBroderick fc Loebll2005^ . 



ACKNOWLEDGMENTS 

We would like to thank Michael Wise for help with the 
scripts for producing FITs images from the calculations. 
H.F. would like to thank the Miller Institute for hosting 
his Visiting Miller Professorship at the Astronomy Depart- 
ment of UC Berkeley, during which significant progress on 
this paper was made. 



REFERENCES 

Aitken D. K., Greaves J., Chrysostomou A., Jenness T., 
Holland W., Hough J. H., Pierce-Price D., Richer J., 2000, 
ApJ, 534, L173 

An T., Goss W. M., Zhao J.-H., Hong X. Y., R oy S., Rao 
A. P., Shen Z.-Q., 2005, ApJ, in press, (|astro-ph/0503527 ) 

Backer D. C., Sramek R. A., 1999, ApJ, 524, 805 

Baganoff F. K., 2003, AAS/High Energy Astrophysics Di- 
vision, 7, 

Baganoff F. K., Bautz M. W., Brandt W. N., Chartas G., 
Feigelson E. D., Garmire G. P., Maeda Y., Morris M., 
Ricker G. R., Townsley L. K., Waher F., 2001, Nature, 
413, 45 

Baganoff F. K., Maeda Y., Morris M., Bautz M. W., Brandt 

W. N., Burrows D. N., 2003, ApJ, 591, 891 
Balick B., Brown R. L., 1974, ApJ, 194, 265 
Belanger G., Goldwurm A., Melia F., Ferrando P., Grosso 
N., Porquet D., Warwick R., Yusef-Zadeh F., 2005, ApJ, 
635, 1095 

Bietenholz M. F., Bartel N., Rupen M. P., 2000, ApJ, 532, 
895 

Blandford R. D., Konigl A., 1979, ApJ, 232, 34 

Bower G. C., Falcke H., Herrnstein R. M., Zhao J., Goss 

W. M., Backer D. C., 2004, Science, 304, 704 
Bower G. C., Falcke H., Mellon R. R., 2002, ApJ, 578, L103 



12 Markoff, Bower & Falcke 



Bower G. C, Falcke H., Wright M. C, Backer D. C, 2005, 
ApJ, 618, L29 

Bower G. C., Goss W. M., Falcke H., Backer D. C., Lithwick 

Y., 2006, ApJ, 648, L127 
Bower G. C., Wright M. C. H., Falcke H., Backer D. C, 

2003, ApJ, 588, 331 
Broderick A. E., Loeb A., 2005, MNRAS, 363, 353 
Brunthaler A., Bower G. C., Falcke H., 2006, A&A, 451, 

845 

Brunthaler A., Bower G. C., Falcke H., Mellon R. R., 2001, 

ApJ, 560, L123 
Corbel S., Nowak M., Fender R. P., Tzioumis A. K., 

Markoff S., 2003, A&A, 400, 1007 
Eckart A., Baganoff P., Morris M., Bautz M. W., Brandt 

W., Garmire G., Genzel R., Ott T., Ricker G., Straub- 

meier C., Viehmann T., Schodel R., 2004, A&A 
Falcke H., 1996, ApJ, 464, L67 
Falcke H., Biormann R L., 1995, A&A, 293, 665 
Falcke H., Biormann R L., 1999, A&A, 342, 49 
Falcke H., Kording E., Markoff S., 2004, A&A, 414, 895 
Falcke H., Markoff S., 2000, A&A, 362, 113 
Falcke H., MeMa F., 1997, ApJ, 479, 740 
Falcke H., Meha F., Agol E., 2000, ApJ, 528, L13 
Fender R., Corbel S. ., Tzioumis T., Mclntyre V., 

Campbell- Wilson D., Nowalc M., Sood R., Hunstead R., 

Harmon A., Durouchoux P., Heindl W., 1999, ApJ, 519, 

L165 

Fender R. P., Gallo E., Jonker P. G., 2003, MNRAS, 343, 
L99 

Gallo E., Fender R. P., Miller-Jones J. C. A., Merloni A., 
Jonker P. G., Heinz S., Maccarone T. J., van der Klis M., 
2006, MNRAS, 370, 1351 

Gallo E., Fender R. P., Pooley G. G., 2003, MNRAS, 344, 
60 

Genzel R., Schodel R., Ott T., Eckart A., Alexander T., 
Lacombe F., Rouan D., Aschenbach B., 2003, Nat., 425, 
934 

Chez A. M., Morris M., Becklin E. E., Tanner A., 
Kremenek T., 2000, Nat., 407, 349 

Chez A. M., Wright S. A., Matthews K., Thompson D., 
Le Mignant D., Tanner A., Hornstein S. D., Morris M., 
Becklin E. E., Soifer B. T., 2004, ApJ, 601, L159 

Ho L. C, 1999, ApJ, 516, 672 

Kording E., Falcke H., Corbel S., 2006, A&A, 456, 439 

Lind K. R., Blandford R. D., 1985, ApJ, 295, 358 

Liu S., Molia F., 2002, ApJ, 566, L77 

Macquart J.-P., Bower G. C, Wright M. C. H., Backer 

D. C, Falcke H., 2006, ApJ, 646, Llll 
Markoff S., Falcke H., Fonder R., 2001, A&A, 372, L25 
Markoff S., Falcke H., Yuan F., Biermanu P. L., 2001, A&A, 

379, L13 

Markoff S., Nowak M., Corbel S., Fender R., Falcke H., 

2003, A&A, 397, 645 
Markoff S., Nowak M. A., Wilms J., 2005, ApJ 
Marrone D. P., Moran J. M., Zhao J. ., Rao R., 2006a, ApJ 
Marrone D. P., Moran J. M., Zhao J.-H., Rao R., 2006b, 

ApJ, 640, 308 
Meha F., 1992, ApJ, 387, L25 
Meha P., Palcke H., 2001, ARA&A, 39, 309 
Merloni A., Heinz S., di Matteo T., 2003, MNRAS, 345, 

1057 

Merloni A., Kording E., Heinz S., Markoff S., Di Matteo 



T., Palcke H., 2006, New Astronomy, 11, 567 
Migliari S., Tomsick J. A., Markoff S., Kalemci E., Bailyn 
C, Buxton M., Corbel S., Fender R. P., Kaaret P., 2007, 
ApJ 

Nagar N. M., Falcke H., Wilson A. S., 2005, A&A, 435, 521 
Nagar N. M., Falcke H., Wilson A. S., Ulvestad J. S., 2002, 
A&A, 392, 53 

Narayan R., Maliadevan R., Grindlay J. E., Popham R. G., 

Gammie C, 1998, ApJ, 492, 554 
Nord M. E., Lazio T. J. W., Kassim N. E., Goss W. M., 

Duric N., 2004, ApJ, 601, L51 
Quataert E., 2002, ApJ, 575, 855 

Reid M. J., Menten K. M., Genzel R., Ott T., Schodel R., 

Eckart A., 2003, ApJ, 587, 208 
Roy S., Pramesh Rao A., 2004, MNRAS, 349, L25 
Schodel R., Ott T., Genzel R., Eckart A., Mouawad N., 

Alexander T., 2003, ApJ, 596, 1015 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Shcn Z.-Q., Lo K. Y., Liang M.-C, Ho P. T. P., Zhao J.-H., 

2005, Nat., 438, 62 
Trippe S., Paumard T., Ott T., Gillessen S., Eisenhauer F., 

Martins F., Genzel R., 2007, MNRAS, 375, 764 
Yuan P., Markoff S., Palcke H., 2002, A&A, 383, 854 
Yuan P., Markoff S., Palcke H., Biermann P. L., 2002, A&A, 

391, 139 

Yuan P., Quataert E., Narayan R., 2003, ApJ, 598, 301 
Yuan P., Shen Z.-Q., Huang L., 2006, ApJ, 642, L45 
Yusef-Zadeh P., Roberts D., Wardle M., Heinke C. O., 
Bower G. C, 2006, ApJ, 650, 189 



